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ABSTRACT 



Context. The Vela Molecular Ridge hosts a number of young embedded star clusters in the same evolutionary stage. 
Aims. The main aim of the present work is testing whether the fraction of members with a circumstellar disk in a sample of clusters 
in the cloud D of the Vela Molecular Ridge, is consistent with relations derived for larger samples of star clusters with an age spread. 
Besides, we want to constrain the age of the young embedded star clusters associated with cloud D. 

Methods. We carried out L (3.78 yum) photometry on images of six young embedded star clusters associated with cloud D of the Vela 
Molecular Ridge, taken with ISAAC at the VLT. These data are complemented with the available HK^ photometry. The 6 clusters are 
roughly of the same size and appear to be in the same evolutionary stage. The fraction of stars with a circumstellar disk was measured 
in each cluster by counting the fraction of sources displaying a NIR excess in colour-colour (HK S L) diagrams. 
Results. The L photometry allowed us to identify the NIR counterparts of the IRAS sources associated with the clusters. The fraction 
of stars with a circumstellar disk appears to be constant within errors for the 6 clusters. There is a hint that this is lower for the most 
massive stars. The age of the clusters is constrained to ~ 1-2 Myr. 

Conclusions. The fraction of stars with a circumstellar disk in the observed sample is consistent with the relations derived from larger 
samples of star clusters and with other age estimates for cloud D. The fraction may be lower for the most massive stars. Our results 
agree with a scenario where all intermediate and low-mass stars form with a disk, whose lifetime is shorter for higher mass stars. 

Key words. Stars: formation - Stars: pre-main sequence - Stars: circumstellar matter - ISM: Vela Molecular Ridge - open clusters 
and associations: general - Infrared: ISM 



1. Introduction 

Open clusters have long been considered as suitable laborato- 
' ries to test stellar evolution theories and sample the Initial Mass 
. | Function (IMF). More recently, Near-Infrared (NIR) cameras 
■ unveiled the precursors of open clusters: young embedded star 
5— i ' clusters. Even more interestingly, it is now clear that a signif- 
5^ , icant fraction of star formation in Giant Molecular Clouds oc- 
curs in young embedded clusters. But less than ~ 10 % of these 
emerge from molecular clouds as bound open clusters (see re- 
view by Lada & Lada 2003). Young embedded clusters are now 
considered as even more reliable sites where to measure an IMF, 
although in this case one has to take into account the bias intro- 
duced by effects such as heavy reddening and the NIR excess 
displayed by young pre-main sequence (pms) stars (see, e. g., 
review by Scalo ll998b . 

Giant molecular clouds hosting a number of young embed- 
ded clusters of the same age represent much more valuable lab- 



* Based on observations collected at the European Southern 
Observatory, Cerro Paranal, Chile, programme 074.C-0630 



oratories, allowing one to test even a higher number of physical 
mechanisms than in single star clusters. In this respect, one of 
the most interesting regions is the Vela Molecular Ridge ( VMR), 
first studied in the CO(1-0) by Murphy & May ( fTWIT ). The 
VMR is a Giant Molecular Cloud complex located in the galactic 
plane, i n the o uter Galaxy (260° < I < 272°, \b\ < 2°). Murphy 
& May d 1 99 It subdivided the region in 4 main clouds, called A, 
B, C and D. The issue of distance was discussed by Liseau et al. 
(fT992b who concluded that A, C and D lie at 700 + 200 pc. The 
molecular gas distribution of the VMR was then studied with a 
higher resolution by Yamaguchi et al. ( 1999), whereas Elia et al. 
( 2007 ) and Massi et al. ( 2007 ) mapped a 1 x 1 sky area towards 
cloud D in the CO(1-0) and 13 CO(2-l) transitions and in the 
continuum 1 .2 mm emission, respectively, with the SEST. The 
dynamical structure of the molecular gas suggests an age of ~ 1 
Myr for this part of cloud D. Star formation in clouds A, C and 
D (both in young embedded clusters and in isolation) was stud- 
ied in a number of works (to quote but a few: Liseau et al. 1992; 
Lore nzetti et al. [T9931 120021 Wouterloo t & Brand [19991 Massi 
et al. [19991 12000 2003; Giannini et al. [20031 [20071 Burkert et 
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al. 120001 Caratti o Garatt i et al [20041 Baba et al.|2 004, 2006; 
Apai et al. l2005l Thi e tal.l2006l De Luc a et al [20071 Ortiz et al. 
120071 Netterfield et al [20091 Olmi et al. |2009l >. 

Massi et al. (2006) exploited the natural "laboratory" pro- 
vided by cloud D by selecting a sample of 6 small young embed- 
ded clusters in a similar evolutionary stage that they observed in 
the NIR JHK S bands. They studied the IMFs and tried to con- 
strain the age, as well. It could be shown that the IMFs are con- 
sistent each others and consistent with a standard IMF (e. g., the 
one proposed by Scalo 1998 ). But the cluster age could only be 
loosely constrained to lie in the interval ~ 1 — 5 Myr. Therefore, 
it is critical to obtain a better age determination. 

A valuable tool to determine the age of young embedded 
clusters is their content of pms stars with a circumstellar disk, 
that are identifiable due to their NIR excess with respect to pho- 
tospheric emission (Haisch et al. 2001 Hillenbrand 2003J. In 
principle, this can be easily achieved by counting cluster mem- 
bers displaying a NIR excess and ones just falling in the red- 
dening band of the main sequence, in a JHKL or HKL colour- 
colour diagram. In the case of the VMR, it is tempting to exploit 
the large number of young embedded clusters of the same age to 
also test whether the fraction of stars with a circumstellar disk 
only depends on the cluster age or other effects have to be taken 
into account. In a scenario where all protostars accrete through 
a disk, the fraction of stars with a circumstellar disk in a cluster 
must depend only on disk lifetimes, hence it is expected to be 
only a function of cluster age, at least in clusters of similar sizes. 
With these two aims in mind, we selected a sample of 6 young 
embedded clusters associated with the part of cloud D mapped 
by Elia et al. (2007 ) and Massi et al. |2007 ), and imaged them in 
the L (3.78 fim) band. The obtained photometry was then com- 
plemented by HK S photometry already available. 

The paper layout is the following: observations and data re- 
duction are described in Sect. [2] the results are discussed in 
Sect.[3]and summarised in Sect.|4] 

2. Observations and data reduction 

We selected the same fields as Massi et al. (2006), excepted 
IRS 19 since the NIR counterpart of the IRAS source is too bright 
in L and would saturate the detector. We replaced it with IRS22, 
another site hosting a young embedded star cluster and associ- 
ated with cloud D. These are the most populated clusters found 
in cloud D of the VMR. The six fields are listed in Tab. [T] along 
with the gas and stellar masses associated with the stellar clus- 
ters. The total gas mass is inferred from Table 3 of Elia et al. 
(2007 ), whereas the total mass in dense gas from Table 1 of 
Massi et al. (2007). We note that the listed total gas masses 
are less than the total masses in dense gas in three cases. The 
gas mass is derived in Elia et al. (2007) by using CO(1-0) and 
13 CO(2-l) observations (and assuming LTE conditions) with 
the highest spatial resolution obtained so far. Nevertheless, their 
13 CO(2-l) observations are spatially undersampled by a factor 
2 and their sensitivity is low (~ 0.7 K), meaning that the com- 
puted masses have to be considered as mere lower limits. In fact, 
Yamaguchi et al. (Q999 ) find much higher gas masses by using 
CO(1-0) and 13 CO(1-0) observations with a larger beam. 

The fields were observed with the camera ISAAC 
(Moorwood et al. [19991 at the ESO-VLT telescope (Cerro 
Paranal, Chile) through the broad-band filter L centred at 3.78 
/mi, in the LW imaging mode. The pixel scale is ~ 0.07 arc- 
sec/pixel and the field of view is 73 x73 arcsec 2 . All observations 
were made in service mode. IRS 17 and IRS 18 were observed 
on 24/01/2005, IRS 20, IRS 21 and IRS 22 were observed on 



25/01/2005 and IRS 16 on 20/02/2005. On each night, images of 
the standard star HD75223 (van der Bliek et al. 119961 ) and dark 
frames were also taken. Sky flats were acquired on 24/01/2005. 

For every field, two sets of 12 chopping sequences of im- 
ages were taken. In each chopping sequence, the pointing was 
switched between two different positions (hereby called on and 
off, although both point towards the target) 30 times, by mov- 
ing the telescope secondary mirror. The integration time in each 
on and off position results from an average of (NDIT=) 9 single 
exposures of (DIT=) 0.11 s. After each chopping sequence, all 
ons were averaged together and the same for all offs, then the 
telescope was nodded, following an ABBAABBA ... sequence. 
The chop throw was 15" and the only difference between the two 
sets of 12 chopping sequences is the throw angle. The first se- 
quence is chopped (in equatorial coordinates) east-west, whereas 
the second one in a direction rotated by 45 degrees with respect 
to the previous. A small jitter (~ 5") was used between ABBA 
chopping sequences. The nodding was in the same direction as 
the chopping, with the same throw. The total integration time per 
field is ~ 24 minutes. 

Data were reduced using IRAfQ routines and following the 
steps outlined in the "ISAAC Data Reduction Guide 1.5", avail- 
able online. Each on and off frame was first dark subtracted and 
then corrected for non-linearity according to Sect. 7.4.2 of the 
Data Reduction Guide (using updated coefficients). After flat- 
fielding, for each single chopping sequence the averaged off 
frame was subtracted from the averaged on frame. All obtained 
on-off frames were then corrected for bad pixels. 

An on-off frame consists of displaced "positive" and "neg- 
ative" images of stars (due to the chop throw). Each AB cycle 
yields four different shifted images of the same field, i. e. two 
on-off and two off-on. For each AB cycle, by subtracting the 
on-off B from the on-off A one obtains a new image where 
residual sky patterns are effectively cancelled. However, there 
are now two symmetric "negative" stars with respect to the cor- 
responding positive ones (twice in counts). This image can be 
multiplied by -1 and shifted to obtain two more independent 
images (also with negative stars). Therefore, we produced 4 sky- 
subtracted, shifted images for each AB cycle. Finally, all sky- 
subtracted images of a same field were registered and averaged 
together, after removal of the "negative" stars, to yield the final 
image. We also constructed a median final image in the same 
way, and an average and a median image for each of the two 
sets with different chopping throw angle. Note that, because of 
chopping and nodding, only the very central area of the final im- 
age exploits the full ~ 24 minutes of integration. The outermost 
image area was exposed for 6 minutes. The effective integration 
time increases from the edge to the centre. 

The seeing was quite good: we found average PSF FWHMs 
of ~ 0.3" (IRS16), ~ 0.5" (IRS17 and IRS18), and ~ 0.4 - 0.6" 
(IRS20, IRS21 and IRS22). Photometry on the L images was 
performed by using DAOPHOT tasks in IRAF. For each field, the 
list of detections obtained with DAOPHOT was visually checked 
on the final median images, to search for both false and missed 
detections. It was very useful having images with different chop- 
ping throw angles to check for stars that could have been can- 
celled in the preliminary on-off subtraction. Aperture photome- 
try was carried out with PHOT, by using an aperture radius ~ 1 
FWHM and inner and outer radii for the sky annulus of ~ 2 and 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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IRAS fl 


RA(2000) 


DEC(2000) 


Total gas mass'' 


Total mass in dense gas'' 


Cluster stellar Mass rf 


source 






(Mq) 


(Mq) 


(Mq) 


IRAS0843 8^4340 (IRS 16) 


08:45:35.8 


-43:51:00 


99 


105 


103 


IRAS08448-4343 (IRS 17) 


08:46:35.0 


^3:54:30 


48 


88 


90 


IRAS08470-4243 (IRS 18) 


08:48:47.7 


^12:54:22 


o 


CO 


89 


IRAS08476-4306 (IRS20) 


08:49:26.6 


^3:17:13 


13 


21 


70 


IRAS08477-4359 (IRS21) 


08:49:32.9 


-44:10:47 


219 


60 


68 


IRAS08485-4419 (IRS22) 


08:50:20.7 


^14:30:41 


o 


(') 


> 38 s 



" also listed the nomenclature used by Liseau et al. (1992); * derived from CO(1-0) and 13 CO(2-l) observations by Elia et al. (2007); c derived 
from continuum 1.2 mm observations by Massi et al. (2007 ); d estimated by Massi et al. |2006); e found CO(1-0) emission (Liseau et al. l 19921 ): ' 
found a dense core at 1.2 mm centred at the cluster (Massi R, private comm.); s inferred from the data in Massi et al. {2003 ). 



~ 4 FWHMs. At last, we performed PSF-fitting photometry with 
ALLSTAR. 

For IRS16, IRS17, IRS 18, IRS20 and IRS21, we comple- 
mented our L photometry with the JHK S photometry (from SofI 
images taken at NTT/ESO) by Massi et al. ( 120061 ). We cross- 
checked both photometries and redid the JHK S photometry on 
few sources that had gone undetected by Massi et al. (2006). In 
the case of IRS 17, the group of sources towards # 40 of Massi et 
al. ([1999), resolved on the ISAAC image, are barely resolved on 
the SofI images. Hence, only for the group of stars close to # 40 
we adopted the JHK S photometry (from ISAAC images taken at 
the VLT) by Giannini et al. (I20051) . 

For IRS22, we performed new HK S photometry on SofI im- 
ages taken on 31/12/2007 as part of ESO programme 080. C- 
0836. These used the Large Field imaging mode (like in Massi 
et al. 120061 ) with a plate scale of ~ 0.29 arcsec/pixel, yielding 
a f.o.v. of 4.9 X 4.9 arcmin 2 . Both at H and at K s , the field was 
imaged in a 2 x 2 square grid, chosen in such a way that all 
frames overlap on a central area ~ 2x2 arcmin 2 containing the 
cluster. For each of the 4 positions, 10 dithered integrations were 
performed, each one a mean of 12 exposures of 2 s. After each 
on-source integration, the telescope was moved to a location off- 
source and a sky image taken with the same averaged exposure 
time. Dome flat field images were taken before the observations. 
Data reduction was carried out by using standard IRAF routines. 
All frames were corrected for cross-talk and flat-fielded. For 
each on-source frame, a sky frame was constructed by median- 
combining the 4 closest (in time) off-source frames after star re- 
moval. All on-source frames were then sky-subtracted and cor- 
rected for bad pixels. At last, they were registered and combined 
together through their average. The L-field is contained in the 
central ~ 2 x 2 arcmin 2 area of the HK S mosaics, so the total 
integration time of interest is 16 minutes. 

Photometry in the HK S band was carried out with 
DAOPHOT tasks in IRAF. We searched for stars with DAOFIND 
and, then, visually inspected the images to remove false detec- 
tions and/or add missed detections (after a cross-check with the 
L image). We used PHOT to perform aperture photometry, with 
an aperture radius of ~ 1 FWHM and inner and outer sky annuli 
of ~ 2 and ~ 4 FWHMs (with median as a sky estimator). The 
PSF FWHM is ~ 0.9". Calibration was obtained by using stan- 
dard stars imaged during the night from the list by Persson et al. 
((19981 . 

Finally, we combined the JHK S and L photometry together. 
Over the total solid angle imaged by ISAAC towards each clus- 
ter, the number of L sources with a K s counterpart found is: 107 
(IRS16), 67 (IRS17), 133 (IRS 18), 49 (IRS20), 71 (IRS21) and 
101 (IRS22). The number of L sources without a K s counter- 
part is instead: 13 (IRS22), 2 (IRS16), 1 (IRS17), 1 (IRS18) 



and (IRS20 and IRS21). The number of K s sources without 
an L counterpart is: 34 (IRS16), 65 (TRS17), 87 (IRS18), 109 
(IRS20), 48 (IRS21) and 105 (IRS22). If K s - L is computed 
for the K s sources without an L counterpart by using the L com- 
pleteness limits estimated in Sect. 13.11 most of these objects fall 
on the right of the main sequence reddening band in a H - K s 
vs. K s - L diagram. I. e., they are consistent with reddened or 
unreddened stars whose L photospheric emission lies below the 
completeness limit (in flux). Some or them might even have a 
NIR excess and still go undetected in L. This also means that the 
available K s photometry is deeper (in terms of sampled stellar 
masses) than our L photometry. 

3. Results 

3. 1 . Completeness limits 

To estimate the completeness limit of our photometry, we first 
constructed histograms of the number of sources vs. L for each 
field, with L binned in 0.5 mag intervals. In all histograms, the 
number of sources increases with increasing magnitude, peaks, 
and then decreases. The peak is the signature of the efficiency in 
finding sources starting to decrease. Using the task ADDSTAR 
in DAOPHOT, we carried out experiments by adding small se- 
quences of artificial stars randomly distributed over each frame. 
We could always retrieve more than 80 % of the artificial stars at 
the central magnitude of each histogram peak. Furthermore, the 
magnitudes measured with aperture photometry (PHOT) were 
within ~ 0. 1 mag of the actual values around the peak magni- 
tude. 

Hence, we estimated completeness magnitudes of L comp i ~ 
14.25 (IRS 16, IRS21 and IRS22), L' ~ 13.75 (IRS 17 and 
IRS 18) and L" ~ 13.25 (IRS 20), that is 1-1.5 magnitudes 
below the detection limits (L ~ 14.5-15.7). The latterones were 
estimated by examining the error vs. magnitude diagrams and lo- 
cating the faintest stars in each diagram. IRS 18 and IRS22, the 
most populous clusters, also exhibit a second peak in the num- 
ber of sources vs. L, located ~ 2 mag below the completeness 
magnitude, that we interpret as an intrinsic feature of the cluster 
population. 

Using the assumed distance modulus to Vela-D (9.225 mag) 
and an estimate of the maximum extinction, we can convert 
^compi C^compi' ^c'ompi-* to a limiting absolute magnitude. In turn, 
this can be converted into a limiting stellar mass if a maximum 
age for the cluster members is also assumed. By using the red- 
dening law by Rieke & Lebofsky ( 1985), we derived complete- 
ness limits in absolute magnitudes of L ~ 4.45 (~ 3.95,3.45) 
for A v ~ 10 mag, and L ~ 2.71 (~ 2.21, 1.71) for A v ~ 40 
mag. To convert these magnitudes into stellar masses, we made 
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use of pms evolutionary tracks by a number of different authors. 
Adopting the tracks by Palla & Stahler (1 19991 1. by a rough com- 
parison with the colours given by Koornneef ( 119831 1 and Bessell 
& Brett (1988 ), the values for Ay ~ 40 mag are consistent with 
M ~ 0.4 Mq (~ 0.6, 1 M Q ) for an age of 1 Myr, and M < 0.1 
M Q (< 0.2, ~ 0.2 Mq) for an age of 0.1 Myr. Adopting the 
tracks by D' Antona & Mazzitelli ( 1994), these numbers become 
M ~ 0.6 Mq (~ 0.8, 1 Mq) for an age of 1 Myr, and M ~ 0.2 
Mq (~ 0.3, 0.4 Mq) for an age of 0. 1 Myr. 

The tracks by Siess et al. (2000), with the colours of Kenyon 
& Hartmann JT995T ), yield M ~ 0.2 Mq (~ 0.3, 0.7 Mq) for an 
age of 1 Myr and M < 0.1 Mq (~ 0.2, 0. 3 Mq) for an age of 0.1 
Myr. At last, the tracks by Baraffe et al. (119981 ) yield M ~ 0.13 
Mq for an age of 2 Myr, if the colours by Kenyon & Hartmann 
(119951 ) are used. However, these colours appear inconsistent with 
those inferred from Baraffe et al. (1998). Using the colours by 
Koornneef (119831 ), one obtains M ~ 0.7 Mq. In summary, the 
L completeness limits correspond to stellar masses in the range 
0.1-0.7 Mq (0.2-0.8, 0.2-1 Mq), seen through an extinction of 
Ay ~ 40 mag. We note that this takes into account the star pho- 
tometric emission only, so lower mass stars are still detectable if 
they exhibit a NIR excess. 

3.2. The observed fields 

Zoom-ins of the inner ~ 1 x 1 arcmin 2 fields imaged in the L- 
band are shown in Figs. |T]to [6] The L images are overlaid with 
contours of the lowest levels of the emission in K s (i. e., start- 
ing from ~ 10cr of the background) measured on SofI (IRS 16, 
IRS 17, IRS 18, IRS20, IRS21; Massi et al. l2006t IRS22, this 
work) images to better show the differences between the two 
bands. Clearly, the diffuse emission in the L band is not as promi- 
nent as in the K s band. The spatial resolution in L also appears 
better than in K s . Furthermore, the brightest K s sources appear 
not saturated in L. All these facts enable us to resolve sources in 
L that are not resolved in K s or lie within intense diffuse emis- 
sion. Note also that we have better spatial resolution, and none 
of the sources is saturated, with respect to Spitzer data. 

A clear example is the central bright source found by Massi 
et al. (120031 ) in IRS 16 and barely resolved in two stars, their # 
90 and # 196. Our L image separates the two objects and further 
shows a fainter star in between them. The two L stars without 
a K s counterpart both lie towards this group of stars, as well. In 
the L image of IRS22, source # 1 1 1 of Massi et al. ( 120031) ap- 
pears actually composed of two close-by stars. Furthermore, four 
of the L stars without a K s counterpart lie towards the compact 
group of sources around #111. 

3.3. The cluster stellar population 

To study the stellar population of the clusters, we computed the 
NIR colours of all sources with L < 14.25 (the estimated com- 
pleteness limit in L for half of the clusters, see Sect. 13. U and 
plotted them in a H - K s vs. K s -L diagram. These diagrams are 
drawn in Fig. [7] for all fields. To discriminate between sources 
with a NIR excess and reddened stars, a knowledge both of the 
colours of unreddened main sequence stars and of the reddening 
law is needed. According to the ISAAC manual, the available L 
filter is more similar to the L' filter described in Bessell & Brett 
(fT98"8l . According to their Table II, K - L and K - L' can dif- 
fer by up to ~ 0.1 mag for main sequence M stars. Conversely, 
the values of H - K are mostly within few ~ 0.01 mag with re- 
spect to those given by Koornneef ( 1983 ). Hence, we adopted the 
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Fig. 1. L-band image of the inner ~ 1 x 1 arcmin 2 region of the 
cluster associated with IRS 16. This is overlaid with contours of 
the emission in the /T s -band (from the SofI image by Massi et al. 
2006). Contours range from ~ 10o" of the background in steps 
of ~ lOcr. The IRAS uncertainty ellipse is also drawn, and the 
most prominent NIR sources identified by Massi et al. (2003) 
are labelled according to the numbers these authors assigned in 
their list. 



colours (and L') by Bessell & Brett ( 119881 ) for the main sequence 
locus drawn in Fig. [7] However, note that the 6 fields were im- 
aged in the K s (SofI) band, that is slightly different from the K 
band adopted by Bessell & Brett (fl9%8l ) or Koornneef (fl983l . 
According to Persson et al. (1998), the scatter in the relation- 
ship between K s - K and stellar colours is mostly due to the 
presence or lack of stellar CO-band absorption that affects the 
K filter more than the K s filter. For the standards listed in their 
Table 2, the differences between K and are at most less than 
few ~ 0.01 mag. Differences between K and K s can be also de- 
rived by using the transformations given by Carpenter (2001) 
for 2MASS. In the colour range spread by main sequence stars, 
the differences between K s (2MASS) and K (in the Koornneef 
or Bessell systems) are always < 0.05 mag. These should be of 
the same order for the SofI filter, so the main sequence colours 
given by Bessell & Brett (1988 ) should be within ~ 0.05 mag 
at most of the actual ones in the SofI system. The adopted red- 
dening law is that of Rieke & Lebofsky ( 119851 ). that Massi et al. 
(2006 ) showed to be consistent with the reddening in the SofI 
JHK S system. 

In all colour-colour diagrams of Fig.|7]it is quite evident that 
part of the sources are aligned within the reddening band of the 
main sequence, whereas a large fraction of sources exhibit a NIR 
excess although they are spread along a direction parallel to the 
reddening vector, as well. The NIR counterparts of the IRAS 
sources (reidentified in Appendix [A} appear usually located at 
the upper edge of this sequence of sources with a NIR excess, or 
of its projection towards the upper-right corner of the diagram 
(e. g., source # 57 in IRS 17 and source # 98 in IRS20), meaning 
they are heavily reddened, as well. In this respect, source # 98 
(IRS20) is by far the most extincted of all counterparts (Ay ~ 
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Fig. 7. H - K s vs. K s - L diagrams for all sources with L < 14.25 (the completeness limit for half of the clusters) in each field. 
Arrows mark lower limits in H - K s . The locus of main sequence stars (Bessell & Brett 1988) is drawn (full line) along with the 
loci of reddened B8 and M5 dwarfs up to Ay - 40 mag (dashed lines) according to the law of Rieke & Lebofsky (119851 I. A few 
remarkable sources discussed in the text are labelled using the numbers in the catalogues of Massi et al. (1999 2003). 



80 mag), consistent with it being viewed through an edge-on 
circumstellar disk. 

There are also datapoints found above the reddening band 
of the main sequence in all fields. We found that most of these 
sources are very close to other sources in the K s frames, so that 



the aperture photometry in is contaminated by the light from 
the nearby stars and the corresponding K s - L values are moved 
leftwards, as well. 

The sources within the reddening band of the main sequence 
exhibit extinctions up to Ay = 20-40 mag. This and the high 
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Fig. 2. L-band image of the inner ~ 1 x 1 arcmin 2 region of the 
cluster associated with IRS 17. This is overlaid with contours of 
the emission in the Zf s -band (from the SofI image by Massi et al. 
2006 ). Contours range from ~ 10cr of the background in steps 
of ~ 10o". The IRAS uncertainty ellipse is also drawn, and the 
most prominent NIR sources identified by Massi et al. (1999 ) 
are labelled according to the numbers these authors assigned in 
their list. 



fraction of sources with a NIR excess, indicate that the fields im- 
aged in L are dominated by the young members of the embedded 
clusters. 

To further prove that the imaged L fields are centred towards 
the cluster cores, we plot the surface density of L sources in 
Fig. [8] for all the fields (left boxes). The surface density is ob- 
tained by counting all sources with L < 14.25 within 14"-side 
squares, displaced by 7" in right ascension and declination from 
each other. The surface density maps are centred on the identi- 
fied NIR counterparts of the IRAS sources (see Appendix [At. 
By binning together all counts per cell towards the 6 clusters, we 
obtained a histogram that is reminiscent of a Poissonian curve 
with 1 average and an excess of counts on the wing. This allows 
us to roughly estimate a standard deviation of 20 stars arcmin -2 
towards all fields. Clearly, all cluster cores fall within the im- 
aged fields. Note that some of the proposed counterparts are not 
located in the centre of the clusters (e. g., source # 57 in IRS 17). 

3.4. The population of NIR stars without an L counterpart 

In Sect. [2] we noted the presence of K s sources without an L 
counterpart in all fields and checked that they are consistent with 
faint stars whose L flux falls below the completeness limit. This 
is confirmed by the magnitude-colour diagrams (K s vs. H - K s ) 
shown in Fig.|9]for all K s sources in the imaged fields without an 
L counterparts. All these sources are actually fainter than K s ~ 
14-15. 

To further confirm their nature, in Fig. [8] we plot contours of 
the surface density of K s sources without an L counterpart (right 
boxes) for all fields, obtained as explained above. Some of the 
fields exhibit a clustering of such sources that is not centred to- 
wards the spatial distribution centre of the L sources. However, 
note that the brightness of most of these K s sources is below the 
completeness limit in K$. Then, part of these faint sources may 
be cluster members, as well, and their different spatial distribu- 
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Fig. 3. L-band image of the inner ~ 1 x 1 arcmin 2 region of the 
cluster associated with IRS 18. This is overlaid with contours of 
the emission in the ^ s -band (from the SofI image by Massi et al. 
2006). Contours range from ~ 10cr of the background in steps 
of ~ lOcr. The IRAS uncertainty ellipse is also drawn, and the 
most prominent NIR sources identified by Massi et al. (1999) 
are labelled according to the numbers these authors assigned in 
their list. 



tion may just reflect variations in the detection efficiency due to 
variations in extinction and diffuse emission over the fields. 



3.5. Contamination of the cluster members by field stars 

To check how much the population of L sources is contaminated 
by field stars, we derived the K Luminosity Function (KLF) of 
each cluster by using all K s sources with an L counterpart. These 
KLFs were then dereddened according to the method discussed 
in Massi et al. (2006 ), and compared with the ones obtained for 
the same clusters by these authors after dereddening and correc- 
tion for field stars (excluded IRS22, not observed by Massi et al. 
2006). Figure [TUl shows that all pairs of KLFs are quite similar 
for all 5 fields at least down to K s = 14, i. e. the completeness 
limit as estimated by Massi et al. (120061 ) for their dereddened 
KLFs. Hence, the L photometry samples the KLFs down to a 
dereddened K s = 14, as well. Some differences are evident at 
the high luminosity ends, where the statistic is however poor. 
A x 2 test on eacn P a ir of KLFs shows that they are equal at a 
significance level > 70 %. This indicates that the L sources fully 
sample the cluster population and that the contamination by field 
stars is very low. 



4. Discussion 

4. 1 . Finding the stars with a circumstellar disk 

The selected clusters are a suitable sample to check whether the 
fraction of members with circumstellar disks depends on only 
the cluster age or other effects have to be accounted for. First, 
a number of facts suggest that all 6 clusters are roughly equally 
old: 
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Fig. 8. Maps of the surface density of stars towards the observed fields. In the left-hand boxes the surface density of sources with 
L < 14.25 is plotted, whereas in the right-hand boxes the surface density of K s sources without an L counterpart is shown. The 
lowest contour (thick line) is 40 stars arcmin~ 2 , the other contours are in steps of 40 stars arcmin~ 2 . Coordinates are offsets from the 
position of the identified counterparts of the IRAS sources. 
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Fig. 9. Ks vs. H - K s diagrams for all K s sources in the observed fields without an L counterpart. The locus of the ZAMS at d — 700 
pc is also drawn, adopting the colours of Koornneef ( 1983 > and the absolute magnitudes given by Allen (119761 I. A few spectral types 
are labelled and a reddening vector Ay = 20 mag is drawn, as well, according to Rieke & Lebofsky d 1 985b . 



1. They are all associated with molecular gas (Massi et al. 1999 
12003] 120061 see Tab. [Tl 

2. they are all associated with dense molecular cores (Massi et 
al. 120071 . see Tab. [2 

3. they are all associated with a FIR point source that has been 
identified (excepted for IRS 16) as an intermediate-mass star 
progenitor (Massi et al. [19991 120031 : 

4. they have roughly the same projected size (Massi et al. 1999, 
120031) : 



5. Massi et al. (2006) showed that the best fits of their KLFs are 
obtained with pms stars of an age ~ 1 — 5 Myr; 

6. Elia et al. (2007) found that four of the clusters are located 
within filaments of molecular gas whose origin is probably 
dynamical with an estimated age of ~ 1 Myr. 

Secondly, all clusters belong to a same region (Vela-D) and 
are, therefore, at a same distance and suffer a similar reddening 
distribution. All of them were observed in JHK S L by using the 
same instruments and instrumental setups. Hence, they are all 
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Fig. 4. L-band image of the inner ~ 1 x 1 arcmin 2 region of the 
cluster associated with IRS20. This is overlaid with contours of 
the emission in the Zf s -band (from the Soft image by Massi et al. 
2006 ). Contours range from ~ 10cr of the background in steps 
of ~ 10o". The IRAS uncertainty ellipse is also drawn, and the 
most prominent NIR sources identified by Massi et al. (1999 ) 
are labelled according to the numbers these authors assigned in 
their list. 

affected by the same environmental and instrumental biases and 
their properties are easier to compare. 

Haisch et al. (120001 Lada et al. (120001 ) and Lada & Lada 
(120031 ) showed that virtually every source with a circumstel- 
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Fig. 5. L-band image of the inner ~ 1 x 1 arcmin 2 region of the 
cluster associated with IRS21. This is overlaid with contours of 
the emission in the /f s -band (from the SofI image by Massi et al. 
2006 ). Contours range from ~ 10cr of the background in steps 
of ~ 10cr. The IRAS uncertainty ellipse is also drawn, and the 
most prominent NIR sources identified by Massi et al. ([1999) 
are labelled according to the numbers these authors assigned in 
their list. 
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Fig. 6. L-band image of the inner ~ 1 x 1 arcmin 2 region of 
the cluster associated with IRS22. This is overlaid with contours 
of the emission in the Zf s -band (from the SofI image taken in 
2008). Contours range from ~ 10cr of the background in steps 
of ~ 5cr. The IRAS uncertainty ellipse is also drawn, and the 
most prominent NIR sources identified by Massi et al. (|2003 ) 
are labelled according to the numbers these authors assigned in 
their list. 
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Fig. 10. Dereddened K luminosity functions for the 6 clus- 
ters obtained by using all K s sources with an L counterpart 
(open triangles, dashed lines), overlaid with the corresponding 
dereddened K s luminosity functions by Massi et al. (2006] full 
squares, solid lines). Errorbars for Poissonian errors are also 
drawn (those for the L sources are slightly shifted along the x- 
axis). The vertical dotted line marks the dereddened complete- 
ness limit (K s ~ 14) estimated by Massi et al. d2006 ) for their K 
luminosity functions. 
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lar disk exhibit a NIR excess in a JHKL diagram. Hence, the 
fraction of young stars with a circumstellar disk can be derived 
by just counting the fraction of objects with a NIR excess in a 
JHKL diagram. We rather used the HK S L diagrams of Fig. |7l 
since many HK S L sources are not detected in J. We adopted the 
main sequence locus of Bessell & Brett i 19881 1 and the redden- 
ing law of Rieke & Lebofsky d 1985b . as discussed in Sect. 13.31 
Conservatively, given an H-K s vs. K s -L diagram we considered 
as sources with a NIR excess those lying below a line parallel to 
the reddening vector, passing through a point 0. 1 mag bluer in 
H - K s than an M5 star. The most optimistic mass limit estimate 
(with Ay ~ 40 mag) corresponding to our completeness L mag- 
nitude, discussed in Sect. 13.11 matches a main sequence M5 star. 
E. g., according to Baraffe & Chabrier ( 1996) the mass of an M5 
V star is ~ 0.13 Mq. In principle, less extincted low mass stars 
(later than M5) could be detected and counted as stars with a 
NIR excess. Adopting a selection cut at a lower L, e. g. at L" v 
reduces by much this kind of contamination. Furthermore, the 
obtained fraction of stars exhibiting a NIR excess is likely to be 
a lower limit. In fact, whereas the H and K s images were usu- 
ally taken under similar seeing conditions, the L ones are char- 
acterised by a better spatial resolution. Hence, K s - L could ap- 
pear bluer than it is for close-by stars, moving the corresponding 
datapoints leftwards in the diagram, i. e., increasing the appar- 
ent number of sources without a NIR excess. Furthermore, all 
sources falling above the reddening band of the main sequence 
have been conservatively counted as stars without a NIR excess. 

As discussed in Sect. 13.51 the contamination by field stars 
is expected to be quite low. It can be estimated by noting that 
our photometry for stars with L < 14.25 appears to sample all 
sources down to a dereddened K^ = 14. By counting all objects 
down to K s = 14 in the dereddened KLFs for field stars esti- 
mated by Massi et al. (2006) and scaling their number to the 
area of the L fields, we derived 5-6 contaminant stars per field. 
This agrees with what found in Sect. 13.51 

The fractions of stars with a NIR excess per cluster are listed 
in Tab. [2 both with and without correction for field stars. The 
quoted uncertainties are obtained by propagating the Poissonian 
errors. We derived the fractions both including all stars down to 
^compi and including only stars down to L" v to test the effect 
of both the different completeness limits and the possible con- 
tamination of the NIR excess region from very low mass stars. 
Clearly, we obtain slightly higher fractions for the lower com- 
pleteness magnitude. If these are also corrected for field star con- 
tamination by using the same numbers as above, they become 
even higher. But the highest values are found for the clusters 
with less sources below £/' ompl (IRS 17 and IRS20) and the con- 
tamination from field stars is now overestimated, as well. Hence, 
the differences for each cluster remain within the errors. 

4.2. Sources of bias in counting stars with a circumstellar 
disk 

There is a major effect that has to be taken into account when de- 
riving the fraction of stars with a NIR excess using cuts in mag- 
nitude. Because of extinction, a cut in magnitude implies that the 
volume in which less massive stars are detected decreases with 
decreasing mass. This would be even worse if the age spread 
were not negligible. If stars of the same mass with and with- 
out a disk were detected within the same volume, even if this 
changed with mass there would be no problems in deriving a 
meaningful ratio. But, because of the NIR excess, stars with a 
circumstellar disk can be detected deeper in the cloud with re- 



spect to stars of the same mass without a circumstellar disk. So, 
the number of stars with a NIR excess is bound to overestimate 
the number of stars with a circumstellar disk used in computing 
the ratio. One can impose a raw limit in mass by selecting only 
sources with a dereddened K s magnitude (see Sect. 13.51 ), A"d ere d, 
less than a limiting value. We computed new fractions by select- 
ing only sources with L < L C ompi and Kd eie d < 14; they are listed 
in Tab. [2] both with and without correction for field star contam- 
ination. The derived values appear to be slightly less than the 
other ones, although within the errors. Because the NIR excess 
is lower in K s than in L, if it is small on average then the degree 
of contamination is much reduced. 

We can take advantage of the fact that we observed similar 
clusters belonging to the same region and with the same instru- 
mental setups, and safely assume that the fraction e of missed 
stars without a disk is the same for all clusters. If «2 is the mea- 
sured number of stars without a NIR excess, the actual number 
of stars without a NIR excess will then be 712/(1 - e). It is easy to 
see that the actual fraction / is given by f m x (1 - e)/(l - f m e), 
were f„, is the measured fraction. Therefore, being /„, roughly 
the same, all derived fractions have to be corrected by a same 
factor. I. e., their comparison is meaningful. And even with a 
pessimistic e = 0.3, /,„ (given in Tab. |2]i should be decreased by 
only 0.85 - 0.9 to obtain /. 

Since the JHK S and L observations have been carried out rel- 
atively far in time, source variability also could affect the derived 
fractions. This is discussed by Lada et al. ((2000), who conclude 
that the fraction of stars falling in the excess region only due to 
variability is at most negligible in the Trapezium cluster. Again, 
we expect that it affects all 6 clusters the same way anyhow. 



4.3. Fraction of stars with a circumstellar disk in the clusters 

Remarkably, 5 out of 6 clusters show the same fraction of mem- 
bers with a NIR excess, within the errors. IRS21 exhibits a lower 
fraction although still within the errors. According to Massi et 
al. ( 120061 . the K s photometry of IRS21 has a K s completeness 
magnitude 1 mag lower than that of the other clusters. A look 
at the colour-colour diagram of IRS21 shows 8 datapoints with 
K s -L > 2 and large photometric errors in H-K s , lying above the 
adopted boundary between stars without and stars with a NIR ex- 
cess. Objects with such a large value of K s -L (i. e., > 2) are usu- 
ally embedded protostars (see, e. g., Haisch et al. 120001 Lada et 
al . 120001 1. All these 8 sources are quite faint in H (with photomet- 
ric errors > 0.2 mag) and undetected in J, so it is likely that their 
H-K s has been overestimated. Dropping these sources from the 
counting increases the fraction of members with a NIR excess 
up to 0.59 (down to L com pi and after correction for field stars), 
whereas counting them as stars with a NIR excess increases the 
fraction to 0.66. In summary, there is evidence that IRS21, as 
well, has the same fraction of stars with circumstellar disks as 
the remaining 5 clusters. 

We also derived the fraction of stars with a NIR excess by 
imposing a lower K^ed limit (i. e., /Tdered < H) in order to se- 
lect the highest mass stars in a cluster. The results are listed in 
Tab. |2] as well. Although the statistic is low and contamination 
by foreground star could be not negligible, the fraction of stars 
with a NIR excess is clearly lower. This would agree with disk 
lifetimes being shorter for higher mass stars. However, spectro- 
scopic observations should be used to infer cluster membership 
for the brightest stars. 
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Table 2. Fraction of HK S L sources with a NIR excess. Also listed the fraction derived by selecting only HK^L sources whose K^, 
after correction for reddening as explained in Sect.|43](/r c | erec j), is less than given thresholds (the selection criteria on L and Adered 
are indicated in the first line). 





L<L" , 

compl 


'-' ^ ^compl 


L < L comp i and A: dcred < 14 


L < Lcompi and A" dcred < 1 1 


Field 


Fraction 


Fraction 


Fraction 


Fraction 


Fraction 


Fraction 


Fraction 


name 


corrected 


uncorrected 


corrected 


uncorrected 


corrected 


uncorrected 


uncorrected 


IRS 16 


0.74 + 0.18 


0.65+0.16 


0.65 + 0.13 


0.60 + 0.11 


0.60 + 0.12 


0.56 + 0.11 


0.46 ± 0.25 


IRS 17 


0.93 ± 0.30 


0.65 ± 0.23 


0.66 + 0.18 


0.57 + 0.15 


0.59 + 0.16 


0.50 + 0.14 


0.56 + 0.31 


IRS 18 


0.73 + 0.13 


0.68 + 0.13 


0.66 + 0.11 


0.63 + 0.10 


0.62 + 0.11 


0.59 + 0.10 


0.43 ± 0.21 


IRS20 


0.91 ±0.28 


0.68 ± 0.23 


0.73 + 0.17 


0.63 + 0.15 


0.66 + 0.17 


0.56 + 0.15 


0.25 ± 0.28 


IRS21 


0.69 ± 0.22 


0.52 + 0.19 


0.49 + 0.12 


0.43 + 0.11 


0.46 + 0.11 


0.41+0.11 


0.29 + 0.16 


IRS22 


0.87 ± 0.24 


0.66 ± 0.20 


0.69 + 0.16 


0.62 + 0.13 


0.67 + 0.15 


0.60 + 0.14 


0.50 ± 0.22 



4.4. Cluster age estimates 

In conclusion, our JHK S L photometry of 6 young embedded 
clusters in Vela-C shows that the clusters host the same fraction 
of stars with a circumstellar disk. This reinforces the notion that 
all stars (at least the low-mass ones) have a circumstellar disk 
at birth that dissipates with time (by accretion onto the central 
star, dynamical interactions, photodissociation and/or evolving 
in a planetary system; for a review of the different mechanisms 
see Hollenbach et al. 12000ft . Interestingly, IRS 16, the only clus- 
ter of the sample associated with an Hn region, is not different 
in this respect. The ionisation source is an early B star (Massi et 
al. 2003) and this suggests that at least not-extreme UV environ- 
ments (i. e., those not originated by O stars) do not significantly 
speed up the dissipation of disks. This agrees with the estimated 
time scales for disk dissipation by UV radiation from external 
sources. The latter is < 10 7 yrs only for the outer region of disks 
(> 10 AU), that does not contribute much to the NIR excess, 
and in Trapezium-like clusters, as shown by Hollenbach et al. 
(120001 

Using the relation found by Haisch et al. (2001), a fraction 
of 0.6-0.7, as observed in the 6 clusters, translates into an age of 
~ 2 Myr, in agreement with what found by Massi et al. (|2006). 
Hillenbrand (2005) derives the fraction of stars with a circum- 
stellar disk in clusters from NIR colours, as well, but in a slightly 
different way. By comparison with her results, we again obtain 
an age for the 6 clusters of 1-2 Myr. Interestingly, the pres- 
ence of an Hn region (263.619-0.533) associated with IRS16 
can be used to constrain the age of the clusters. According to 
Caswell & Haynes (1987|, the size of the radio emission at 5 
GHz is 1'. This agrees with the size measured in the K s im- 
age, that is ~ 2'. Assuming a radius of 1', this corresponds to 
0.2 pc at a distance of 700 pc. Massi et al. (120031 1 found that 
the ionising source is probably a B0-B2 V star. The Stromgren 
radius of an Hn region originated by a B0-B1 ZAMS star is 

—2/3 

r c = 5.17 - 0.75 x n 3 pc, where «3 is the electron density 
in units of 10 3 ctrT 3 (Church well & Walmsley 119731 1. Then, 
263.619-0.533 appears not to have reached the Stromgren ra- 
dius, unless it expanded in a dense (10 4 -10 5 cm 4 ) gas now 
cleared away. Even in the case of very dense gas (10 5 cm -3 ), the 
Hn region would be greater than the Stromgren sphere only for 
an ionising star later than B 1 V. Before reaching the Stromgren 
radius, an Hn region expands at a speed exceeding the sound 
speed in an ionised medium, i. e. > 10 km s (see e. g. Dyson 
& Williams [T997j l. Hence, 263.619-0.533 would have expanded 
to its current radius in < 2 x 10 6 yrs. Even in the case of a Bl 
star and a dense (10 5 crrT 3 ) medium, the dynamical time of the 
Hn region would be short (in comparison to the time needed to 
reach the Stromgren radius), i. e. 4 x 10 4 yrs (Dyson & Williams 



I1997K Then, the dynamical age of the Hn region is anyway fully 
consistent with the inferred cluster age. 

Given that the age of the clusters is constrained at ~ 10 6 yr, 
it would be noteworthy to find out which is the physical process 
that can disperse the circumstellar disks around the sources in 
such a time scale. As for low mass stars, of the mechanisms dis- 
cussed by Hollenbach et al. (2000) only viscous accretion onto 
the central stars seems to be efficient in dispersing the inner disks 
so rapidly. Stellar encounters are irrelevant for small clusters, the 
time scale being at least ~ 10 8 yrs according to Hollenbach et al. 
(2000) for the star densities measured by Massi et al. (f2006). 
This is confirmed, e. g., by the simulations of Adams et al. 
(120061 

5. Conclusions 

We have carried out imaging in the L band (3.78 /im) with 
ISAAC at the VLT, of 6 young embedded star clusters associ- 
ated with cloud D of the Vela Molecular Ridge. These are in a 
similar evolutionary stage and, therefore, of the same age. The 
L photometry is complemented with HK S photometry from SofI 
images. We used H-K s vs. K s -L diagrams to derive the fraction 
of cluster members with a circumstellar disk. Our L images have 
better spatial resolution and no saturation problems in compari- 
son to Spitzer/IRAC images. The main sources of bias and errors 
are carefully discussed. The main points of the present study are 
summarised as follows: 

1 . We revised the NIR infrared sources that had been associated 
with the IRAS point sources in each field. 

2. By selecting K s sources with L counterparts, we constructed 
the cluster KLFs. Once dereddened by following the method 
outlined by Massi et al. (2006), we showed that these are 
statistically consistent with those obtained by Massi et al. 
(2006). This demonstrates that the L images provide signif- 
icant samples of the cluster populations and allowed us to 
estimate the degree of contamination from field stars. 

3. Using different selection criteria (also based on the dered- 
dened KLFs), we estimated a fraction of stars with a circum- 
stellar disks ~ 0.6 - 0.7 that appears to be the same for all 
clusters, within errors. This shows that the clusters host the 
same fraction of stars with a circumstellar disk. I. e., the lat- 
ter depends only on cluster age. 

4. We showed that the impact of the main sources of bias is 
much reduced when observing clusters belonging to a same 
molecular cloud with the same instruments and instrumental 
setups, as is the present case. 

5. The derived fraction of stars with a circumstellar disk con- 
strain the age of the 6 clusters to 1 - 2 Myr, refining previous 
age estimates. 
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6. The fraction of stars with a circumstellar disk appears to 
be lower for the most massive stars. This would agree with 
shorter disk lifetimes in massive stars, although the results 
could be biased by a poor statistic. 
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Appendix A: Individual sources 

The observed fields were originally selected by Liseau et al. 
i 1992| l, as sites of IRAS point sources with red colours. Massi et 
al. ([T999 2003) found that the fields are associated with young 
embedded clusters and that all IRAS uncertainty ellipses lie to- 
wards the centres of each cluster. These authors suggest that 
most of the FIR luminosity from each IRAS point source is con- 
tributed by protostars whose NIR counterparts lie within (or very 
near to) the uncertainty ellipse and become brighter and brighter 
with increasing wavelengths. Hence, they should become quite 
evident in the L band. To check this, the NIR SEDs of the 
brightest L sources in each imaged field are plotted in Fig. IA.ll 
Clearly, all plots (one per field) exhibit at least one source with 
a steep SED whose flux increases with increasing wavelengths 
and dominates over the other objects in L. The HK S L magnitudes 
were converted into fluxes following Megessier (1995). A brief 
description of the 6 fields follows. Each field will be identified 
by using the IRS notation of Liseau et al. £1992). 

IRS 16 is different from the other 5 IRAS sources in the 
present sample, in that it was not considered as a 'bona fide' 
Class I source by Liseau et al. (1 19921 1. In fact, it is the only one 
associated with an Hn region, i.e. 263.619-0.533 (Caswell & 
Haynes I1987I I. Massi et al. (120031 ) show that the IRAS uncer- 
tainty ellipse is roughly centred on a B0-B2 star (their # 90, see 
also Fig.[TJ. Massi et al. (2007) found that the IRAS uncertainty 
ellipse lies towards a minimum in 1 .2-mm continuum emission, 
in the middle of three prominent peaks (their MMS1, MMS2 
and MMS3). It is possible that the FIR emission is originated by 
a UCHn region close to the central B star. The source exhibit- 
ing the steepest rising SED in Fig. IA.1I is # 123 of Massi et al. 
(2003|, but it lies outside the uncertainty ellipse. Another NIR 
source with a rising SED is # 100 of Massi et al. (2003), that lies 
within the IRAS ellipse (see Fig.[TJ. Both were already identified 
in JHK by Massi et al. d2003l . 

The NIR counterpart of IRS 17 is clearly the one already 
proposed by Massi et al. (1999), i. e. their # 57. Giannini et 
al. d2005) showed that this NIR source lies towards a dense 
molecular core found in 1 .2-mm continuum emission (MMS4 of 
Massi et al. 2007). They also searched for the driving source of 
a prominent jet, that lies towards a group of NIR sources outside 
the IRAS uncertainty ellipse, on the west. These coincide with 
source # 40 of Massi et al. ([T999 see also Fig. [2), resolved into 
more than one star on ISAAC NIR images presented in Giannini 
et al. (120051 1. Of these, # 40-2 is bright at L and exhibits a rising 
SEDs, as well (Fig. [01 

As for IRS 18, our L image again confirms that the NIR coun- 
terpart of the IRAS source is the one already identified by Massi 
et al. d 19991 1. i. e. their # 1 19. This exhibits a risi ng SE D and is 
one of the brightest L sources in the field (see Fig. lA.U . 

No doubts that the NIR counterpart of IRS20 is the one al- 
ready identified by Massi et al. (1 19991 1, i. e. their # 98. It displays 
a steeply rising SED and is by far the brightest L source in the 
field (see Fig. lA.U . Figure|4]evidences that # 98 lies at the centre 
of a bipolar nebula. Giannini et al. d20071 > found a bipolar jet in 
the H2 2.12-yum emission, centred at # 98 and aligned with the 
nebula, suggesting that it is driven by # 98. This is also consis- 
tent with this source being viewed through a disk edge-on, as 
suggested by Massi et al. (2003]!. 

Instead, Fig. IA. II rules out that # 50 of Massi et al. ( 1999) 
is the NIR counterpart of IRS21, as proposed by these authors, 
lying within the IRAS uncertainty ellipse. The most suited can- 
didate appears to be # 27 of Massi et al. (1999), that however 
is located outside the IRAS ellipse (see Fig. |5j. Source # 35 of 



Massi et al. ( 1999) exhibits a rising SED, as well, but lies outside 
the IRAS ellipse, as well. 

Figure lATI also suggests that the NIR counterpart of IRS22 
might not be # 1 1 1 of Massi et al. (2003 ), but its photometry suf- 
fers from the fact that # 1 1 1 is actually composed of two stars, 
unresolved in H and K s . The most likely NIR counterpart is in- 
stead # 115 of Massi et al. (|2003), that exhibits a steeply rising 
SED and falls within the IRAS uncertainty ellipse (see Fig. [3J. 
However, we cannot rule out that one of the two sources com- 
posing #111 has a steeply rising SED. 
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Fig. A.l. Spectral energy distributions for the brightest L sources found in all 6 fields. The field name is on top of each box. The 
steeply rising SEDs are marked by thick lines. The most remarkable sources are labelled according to the catalogue numbers by 
Massi et al. (fT999l l2003t. 



